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he second-order 
"^^^^^ system is common 

in nearly all branches 
of engineering — car sus- 
pensions, robotic servos, audio filters, 
to name a few — and so, a good under- 
standing of it can be a powerful tool 
in your arsenal. But, the approach 
taken in many engineering schools 
discourages students and leaves them 
with an understanding that associates 
black magic with the topic. 

In this series, I'll dispel that notion 
by demonstrating how common sec- 
ond-order systems are and by showing 
how analogies can be drawn among 
them. Yes, I'll get into mathematics, 
even calculus, but I promise that it will 
be painless. And when you notice how 
widely applicable the math is, you 
won't mind at all. 

In Part 1, I'm going to talk about the 
basic second-order system in terms of a 
physical model (mass-spring-damper| 
and an electrical model (LCR circuit). 
You'll learn the equivalencies between 
these two models and derive the basic 
mathematics that describe them. 

Part 2 will expand the discussion to 
the use of second-order systems as 
filters and then move on to cover 
servomechanisms. You'll learn how 
negative feedback affects the behavior 
of the system and extends the math- 
ematical model. Other implementa- 



tions based on op-amps and DSP (digi- 
tal signal processing) will be dis- 
cussed, too. 

Part 3 will bring the discussion 
home with a detailed analysis of some 
real-world servomechanisms such as a 
mechanical positioner for a robot and 
an electronic phase-locked loop. I'll 
discuss how to predict their perfor- 
mance. How unavoidable non- 
linearities affect the behavior of the 
system and ways to mitigate their 
effects will be covered. 

WEIGHTS, SPRINGS, AND DAMPERS 

In general terms, a second-order 
system is any system that has a mea- 
surement that can change. It stores 
energy related to that change (restor- 
ing force), stores energy related to the 
rate of change (inertia), and dissipates 
energy and comes to a stop (resis- 
tance). Nearly every textbook about 
this subject starts by discussing the 
physical system of mass, spring, and 
damper, so I'll start there. 

Figure 1 shows the arrangement; a 
movable mass is attached to both a 
spring and dashpot (shock absorberl, 
and the other ends are fixed. The verti- 
cal position of the mass is the change- 
able measurement. The spring provides 
the restoring force and stores potential 
energy in proportion to how far the 
mass is moved from its equilibrium 
position. The mass stores kinetic en- 
ergy in proportion to how fast it is 
moving. The dashpot dissipates energy 
(turns it into random heat) at a rate 
proportional to the velocity of the free 
end. Assume that all of the friction in 
the system is embodied in the dashpot. 

The mass's important parameter is, 
well, its mass in units of kilograms. 
(I'll stick to the metric system 
throughout this series because it re- 
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Figure 1— The basic second-order mechanical 
system consists of a spring, mass, and damper, or 
dashpot. 



Under- 
standing 
second- 
order systems isn't a 
mystery, but with so 
many options and ap- 
plications, it's hard to 
get an understanding 
of the principles. So, 
Dave kicks off this 
series by explaining 
some of the basics. 
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Figure 2— A system of just a mass and a damper 
quickly reaches a terminal velocity, where the accelera- 
tion has decayed to zero. 

quires fewer conversion constants.) In 
the earth's gravity field, which has an 
acceleration of roughly -10 m/s 2 , it 
takes a force of 10 newtons (N) to 
suspend a mass of 1 kg in midair. 

The spring is characterized by its 
spring constant (KJ, which indicates 
how much restoring force it develops 
for a given amount of displacement. 
The metric units would be newtons 
per meter. The restoring force is in the 
opposite direction relative to the dis- 
placement; hence, if you move the 
mass downward, the spring's force is 
upward, and vice versa. 

Assume the spring is -100 N/m. If 
you hang your 1-kg mass on it, the 
spring will need to develop a 10-N 
upward force to counteract the -10-N 
(downward) force of gravity. As a 
result, it will stretch by: 



displacement - - 10 N .-0.1m- 
K s -100 N/m 



-10 cm 



The dashpot also develops a nega- 
tive force, but this force is propor- 
tional to the velocity of the free end 
rather than its position. The dashpot is 
characterized by the constant K d in 
newtons per meter per second. The 
definition of "work" is a force operat- 
ing over a distance. It has the same 
units as energy, which are joules (J), and 
1 J is defined as a force of 1 N operat- 
ing over a distance of 1 meter, or 1 
newton-meter. 

The work associated with moving 
the free end of the dashpot against its 
force is turned into randomized heat 
and is essentially lost to the system. 
This is in contrast to the spring, 



which stores that work as potential 
energy and returns it to the system 
when the movement reverses direc- 
tion. Let's play with different values 
of K d and determine what effect this 
has on the response of the system. 

Imagine the dashpot filled with 
thick, syrupy oil. It requires a great 
force to move it even slowly; it has a 
high K d , say, 10 N/m/s. If you hang 
the 1-kg mass from the dashpot, it 
quickly reaches a terminal velocity of 
-0.1 m/s. The work, or energy, needed 
to move the dashpot comes from the 
fact that the mass is descending in 
earth's gravitational field. Terminal 
velocity occurs when the two forces 
balance; the 10 N of gravitational 
force on the mass is balanced by 10 N 
of upward force from the dashpot 
when the velocity is -0.1 m/s. 

A SIMPLE SIMULATION 

Let's put together a mathematical 
simulation. I use MathCad because it 
makes it easy to see the math and 
graphics at the same time, but you 
can use a spreadsheet (see "Simulat- 
ing in a Spreadsheet" sidebar) or your 
favorite programming language. 

The general technique is simple: 
set up three variables that represent 
position, velocity, and acceleration, 
set their initial values, and write 
simple equations that update the 
three values after a short increment of 
time. After a few hundred (or a few 
thousand) timesteps, you can graph 
the evolution of the system by graph- 
ing the three values. In MathCad, use 
a vector to hold the three values so 
that you write a single equation to do 
the updates for each timestep. 

Start by entering the system 
parameters (gravitational accel- 
eration, mass (M), and damper 
parameter, respectively): 



G--i0x m 
M= 1 xkg 

K d — 10 x- N 



(¥) 

Also, set your timestep and num- 
ber of steps and set up an index 
variable. The following equations 
are timestep, number of steps, 
and index variable, respectively: 



dT = 0.001 x s 
N = 200 
1= l...n 

The next step is to enter the initial 
conditions (position, velocity, accel- 
eration, respectively): 

p - 0.1 x m 
A = G 

Next, enter the expressions that 
control updating. Position is easy, it 
gets incremented by the previous ve- 
locity value multiplied by the 
timestep. Similarly, velocity gets 
incremented by the previous accelera- 
tion value multiplied by the timestep. 
Acceleration is trickier, it's the result 
of the net sum of the forces (in this 
case, the constant acceleration of grav- 
ity and the damper force that is pro- 
portional to the velocity) divided by 
the mass. Use vector notation to en- 
sure that all three values get updated 
at each timestep before moving on to 
the next, putting the position on top, 
velocity in the middle, and accelera- 
tion (a) on the bottom: 



Pi -i +4 Txv i-' 
Vj_, +dTxa,_, 

M 



Figure 2 shows the results, with 
the mass quickly reaching a constant 
terminal velocity when the force of 
the dashpot balances that of gravity. 
You can also think of this in terms of 
energy flow,- terminal velocity occurs 
when the work being put into the 
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Figure 3— An overdamped system responds to distur- 
bances with a slowly decaying exponential movement. 
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Current 


Amperes (coulombs/second) 


Physical force 


Newtons 


Electromotive force 


Volts 


Physical work 


Newton-meters 


Electrical work 


Volt-coulombs 


Mass 


Kilograms 


Coil 


Henries 


Spring 


Newtons/meter 


Capacitor 


1 /farads (volts/coulomb) 


Dashpot 


Newtons/meter/second 


Resistor 


Ohms (volts/coulomb/second) 



Table 1— There is a direct analogy for each parameter of the physical system in the electrical system. 



system by the gravity field equals the 
energy being dissipated by the dashpot 
as heat. Remember this point, because 
I'll come back to it later. 

Now, let's go back to the system 
that includes the spring. It has an 
equilibrium point defined by where 
the weight of the mass is balanced by 
the upward pull of the stretched 
spring. If you lift up the mass from 
this point and then release it, it will 
gradually return to that equilibrium 
point. The acceleration of the mass is 
no longer constant, because as it 
moves down, the spring exerts an 
increasing upward force that partially 
cancels the force caused by gravity. 



The other upward force on the 
mass comes from the resistance of the 
dashpot. The work being done on the 
mass by the spring (and gravity) must 
equal the work being done on the 
dashpot by the mass and getting dissi- 
pated. But, the dashpot's force is pro- 
portional to the velocity of the mass 
and the spring's force is proportional 
to the position of the mass. The net 
force on the mass is proportional to 
its acceleration. 

So, the position, velocity (time 
derivative of the position), and accel- 
eration (time derivative of the veloc- 
ity) of the mass decay to zero in 
proportion to each other. One curve 



whose derivative var- 
ies in proportion to 
itself is the exponen- 
tial (the general shape 
of the path followed 
by the mass). 

To see the effect of 
adding the spring, add 
its parameter at the 
top of the simulation 
you did previously and 
modify the updated equation to calcu- 
late the spring parameter: 



-100 x 



N 



_ , + dT x 
_ , + dT x a, 
x v. - 1 + K. 



M 



You'll need more steps in order to see 
everything else, so change N to 5000 
(see Figure 3). 

RESONANCE 

Now imagine a dashpot that's 
filled with only air with K d of 2. If you 
change the parameter in the simula- 
tion and rerun it, you get the behavior 
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shown in Figure 4. Because the dash- 
pot provides so little force, the mass 
has quite a bit of momentum built up 
when it reaches the equilibrium posi- 
tion, so it keeps going until the in- 
creasing force from the spring 
becomes enough to stop it. 

However, now it's too low, so it 
starts back in the other direction. 
This process keeps repeating at a 
regular rate, called the natural, or 
resonant, frequency of the system. A 
system configured this way is said to 
be underdamped. The oscillations do 
die out eventually, because- the dash- 
pot still is dissipating energy from 
the system. 

The resonant frequency (F ) is 
primarily determined by the mass 
and spring (although, as I'll explain 
later, the dashpot does affect it). A 
higher K s will raise the frequency and 
a larger mass will lower it. Ignoring K d 
for the moment, the resonant fre- 
quency is: 

With the values just mentioned, this 
equals 1.592 Hz. In fact, you can see 
about eight cycles in the 5-s elapsed 
time of the simulation. 



Time 
(seconds) 
0.0 1.0 2.0 3.0 4.0 5.0 




Figure 4— An underdamped system oscillates in response to 
disturbances at a particular resonant frequency. This oscilla- 
tion dies out exponentially. 



You can imagine a case somewhere 
between these two extremes of K,. At 
some point, as the resistance of the 
dashpot is gradually decreasing, the 
system's response is just about to 
become periodic. In other words, 
there's a definite boundary between 
the overdamped and underdamped 
regimes. When the system is in this 
configuration, it is said to be critically 
damped. I'll leave this subject for 
now, but I'll return to it later. 



COILS, CAPACITORS, AND RESIS- 
TORS 

Now let's take a look at the ca- 
nonical electrical second-order sys- 
tem, which consists of a coil, 
capacitor, and resistor (see Figure 5). 
The coil is equivalent to the mass, the 
capacitor is equivalent to the spring, 
and the resistor corresponds to the 
dashpot. In fact, for every parameter 
or quantity in the physical system, 
there's a corresponding one in the 
electrical system (see Table 1). 

Start by defining the electrical 
equivalent of position of the mass to 
be charge (on the capacitor), which is 
measured in coulombs. This means 
that the electrical equivalent of veloc- 
ity, or meters per second, is coulombs 
per second, or amperes (A). This 
seems natural enough. 

Next, define the electrical equiva- 
lent of mechanical force, or newtons, 
to be electromotive force, or volts. 
Using both of these definitions, the 
mechanical concept of work = force x 
distance leads to the electrical con- 
cept of work = volts x coulombs. Fol- 
lowing this to its natural conclusion, 
power is defined as work per unit of 
time. If you write power - volts x 
coulombs/seconds and note that cou- 



Simulating in a Spreadsheet 

Even if you don't have a mathematical modeling soft- 
ware package like MathCad, you can do this type of simu- 
lation with a spreadsheet program. I'm going to briefly 
outline the general technique, which should work with 
nearly any spreadsheet program. 

The idea is to dedicate one column of the spreadsheet 
to each variable in the simulation. In the top row, label a 
column for each variable. Put the initial value of each 
variable in the next row, below its label. It's a good idea 
to have a column that represents time, with an initial 
value of zero. The other columns will represent accelera- 
tion, velocity, position, and so on. 

After you have written the mathematical expressions 
that update the variables, enter them as formulae into the 
third row of the spreadsheet and use the values of the 
system variables from the row above. The update formula 
for time is simply the previous value plus a constant that 
represents the time interval for each step of the simula- 
tion. 

Here's the trick: using whatever mechanism your 
spreadsheet uses, select the cells in the third row and 
replicate them downward, using as many rows as the num- 
ber of timesteps you want to simulate. Make sure that the 



cell references for system variables remain relative, 
always using values from the row immediately above. 

Now you should see a huge field of numbers that show 
the evolution of the system over time. If your spreadsheet 
supports the generation of graphs (nearly all of them do), 
you can graph the columns representing the system vari- 
ables against the column representing time, and you 
should get graphs similar to the ones that MathCad gener- 
ates. 

If you want to change the initial conditions, modify 
the values in the second row of the spreadsheet. If you 
want to change the update formulas, do that in the third 
row, but then you'll need to replicate the cells below the 
one(s) you change again. 

You can avoid this second replication if you know how 
to generate absolute as well as relative cell references in 
your formulae. Instead of putting constants into the for- 
mulae, replace them with absolute references to cells near 
the top of the spreadsheet. Then if you want to change the 
parameters of the simulation, you can just change the 
contents of those cells and the rest of the spreadsheet will 
update itself. 
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Figure 5— The basic second-order electrical system 
consists of a capacitor, coil, and resistor. 

lombs/seconds = amperes, you get 
power = volts x amperes, which we're 
all familiar with. 

The remaining equivalencies among 
the mechanical components and their 
electrical counterparts may look 
strange at first, but let's work through 
the details. 

The coil, whose inductance is mea- 
sured in henries |H), is equivalent to 
the mass, which is measured in kilo- 
grams. Using the relationship F = M x 
a, kilograms can be thought of as F/a, 
or newtons per meter per second 
squared. By definition, a coil is a device 
that develops a voltage across itself in 
proportion to the rate at which the 
current is changing, which means that 
henries can be thought of as volts per 
ampere per second, or volts per cou- 
lomb per second squared. You end up 
with the corresponding fundamental 
units you started with in correspond- 
ing positions in the two expressions. 

The electrical capacitor is equiva- 
lent to the mechanical spring. You 
measured K s in units of newtons per 
meter. In a capacitor, Q (charge) = C 
(capacitance) x V (voltage), so C = Q/V. 
In order to establish the equivalence 
between K s and the capacitor, you need 
to use 1/C, or units of 1/farads. In ef- 
fect, a stiff spring (large KJ corresponds 
to a tiny capacitor (low capacitance). 

Finally, the electrical resistor corre- 
sponds to the mechanical dashpot. 
Using the relationship V = I x R, ohms 
can be thought of as V/I, or volts per 
coulomb per second. Again, the units 
match up in the expressions for the 
two systems. 

Everything I said about the me- 
chanical system also applies to the 
electrical system when you substitute 
the correct quantities and units. The 
first simulation you did corresponds to 
a circuit with just a coil and resistor. 
The force of gravity is equivalent to 



putting a battery in series, therefore, 
you can see how the terminal velocity 
of the mechanical system corresponds 
to the VDC current that flows in the 
R-L circuit. 

For the second simulation, you put 
the capacitor back in the circuit. It 
now limits the flow of charge in the 
same way that the spring limits 
movement. And in the third simula- 
tion, lowering K d is equivalent to 
lowering R, and then the electrical 
circuit rings. When you change the 
units in the formula for the resonant 
frequency, it becomes: 



F ° = 2SV (lITc) " 2k 



xvl x C 



SHAKE, RATTLE, AND ROLL 

The next question is, what happens 
if you apply a continuous stimulus to 
either the physical system or electrical 
system, instead of just nudging it and 
watching what happens? For a number 
of reasons, the sinewave is the sim- 
plest form of continuous motion, so 
let's study that. Fourier states that you 
can make any other periodic motion 
by adding sinewaves, which makes it' 
a good place to start. 

Start with the physical system. You 
can apply a stimulus by making the 
fixed anchor movable and then either 
applying a force to it or moving it a 
specific amount. In either case, moving 
the anchor requires that an entity 
outside of the system do work — an 
energy input to the system. The only 
energy output from the system still 
occurs via the damper (or resistor). 
Therefore, it seems obvious that the 
system response will move toward a 
condition in which the dissipation of 
the damper equals the input energy. 

Clearly, both the spring and dash- 
pot will attempt to transfer any mo- 
tion of the anchor to the mass, whose 
inertia will resist it, causing the 
spring and dashpot to change length. 
If you do this slowly, the spring and 
dashpot don't change length by much; 
instead, the entire system, including 
the mass, simply moves up and down 
in sync with the anchor. 

On the other hand, if you do it 
quickly, the inertia of the mass holds it 
nearly still and all of the motion is 
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Figure 6— The system response under continuous stimulation 
shows both a transient part and steady-state part. 



taken by the spring and dashpot. Be- 
tween these two extremes is where 
more interesting things happen. 

Applying a force to the anchor is 
equivalent to a series voltage source. 
Actually moving the anchor is equiva- 
lent to using a current source in paral- 
lel with the inductor. The analogy is 
that if the mass (coil) is large, it moves 
(carries current) hardly at all, while all 
of the motion (current) is taken by the 
spring (capacitor) and damper (resistor). 

If the frequency of the voltage or 
current source is low, the voltage (and 
charge) on the capacitor will track that 
of the source closely. But if the fre- 
quency is high, it will change little. 

In an overdamped system with a 
large K J( most of the force applied to 
the anchor will be coupled with the 
mass via the damper rather than the 
spring. The electrical analogy is that 
the high R, rather than the impedance 
of the capacitance, limits the move- 
ment of charge through the inductor 
and the rest of the circuit. 

In the underdamped system, differ- 
ent things happen depending on 
whether the applied frequency is less 
than, equal to, or greater than the 
system's resonant frequency. Let's 
update the simulation to incorporate 
this external force. 

An external force applied to the 
anchor point gets added to the other 
forces acting on the mass. Set the 
initial position equal to the equilib- 
rium position to simplify the interpre- 
tation of the results (see Figure 6). 

The first three equations are fre- 
quency, anchor force, and initial posi- 
tion, respectively: 



You can see irregularities in 
the graphs near the beginning, 
before things settle down to a 
regular pattern. Looking at the 
acceleration graph, you can see 
that the transient part has a 
frequency component that is 
similar to the natural frequency 
of the system calculated earlier and an 
shape that resembles the graphs in 
Figure 4. This resonance is excited by 
the fact that the simulation applies 
the acceleration suddenly at time zero 
(a phenomenon known as jerk). After 
this transient response dies, the only 
frequency component present is the 
0.5 Hz of the stimulation. 

UNTIL NEXT TIME 

You've observed the analysis of 
both the transient and steady state 
responses of some simple second-order 
systems. It's been an interesting men- 
tal exercise, but I'm sure you're won- 
dering how all of this can be applied 
in the real world. Well, don't worry, 
that's what I'll get into in Part 2. S 

Dave Tweed is an independent con- 
sultant. He has been developing hard- 
ware and real-time software for 
microprocessors for many years. His 
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computer design from supercomput- 
ers to workstations, digital telecom- 
munications systems, and the appli- 
cation of embedded microcomputers 
and DSPs. You may reach him at 
dtweed@acm.org. 
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David Tweed 



Fundamentals of 
Second-Order Systems 

Part 2: The Tools of the Trade 



W 



elcome to Part 
2. I'm presenting 
some of the common 
theoretical background 
and mathematical tools for many of 
the dynamic systems found in all 
branches of engineering. These include 
filters of various types and servo- 
mechanisms such as phase-locked 
loops and robotic positioners. 

Part 1 covered the basic setup of 
both mechanical and electrical second- 
order systems. The former consists of a 
mass, spring, and damper, and the lat- 
ter consists of a coil, capacitor, and 
resistor. Every quantity in each system 
has an analogue in the other, and the 
same math is used to describe both 
systems. A simple numerical simula- 
tion was also set up. 

Now it's time to take a deeper look 
at the underlying mathematics and 
derive more powerful ways of describ- 
ing such systems. We'll take a look at 
how these systems can be used as fil- 
ters, and then what happens when 
negative feedback is added to the sys- 
tem, as in servomechanisms. 

I'm not going to spend a lot of time 
talking about the Mathcad simula- 
tions, but there are files that you can 
download and play with for most of 
the figures in this installment. 



MATHEMATICAL MODEL 

The mathematical formulae you 
built so far are just stepwise numerical 
simulations. But, what about a closed- 
form, continuous solution? It turns out 
not to be too difficult, and involves 
only a little calculus. Bear with me and 
we'll get through it quickly. 

Start by listing the constraints of 
the components and constraints arising 
from how they're connected. The coil 
develops a voltage in proportion to its 
inductance and the time derivative of 
the current through it: 

V IO „=LdI 

The capacitor carries a current that is 
proportional to its capacitance and the 
time derivative of the voltage across it: 



dV„„ 
I = C-j^ 
dt 



Or equivalently, the capacitor voltage 
is the integral of the current through it 
divided by the capacitance: 



-if 



Ixdt 



The resistor develops a voltage in 
proportion to its resistance and the 
current through it: 

V = R x I 

res 

The current, /, is the same every- 
where in the circuit. 

The voltages must sum to zero 
around the circuit: 



V + V 

cap n 



+ v oll -o 



Or, after making these substitutions: 



i! 



Ixdl +RxI + L^ =0 
at 



This final equation is the basic second- 
order problem. If both sides of it are 
differentiated once with respect to 
time, you get a linear, second-order 
homogeneous, differential equation: 

1. + r£ + l£i=o 

Solving such equations is difficult. 
However, note that here you need a 
function that, when added to its first 
and second derivatives, sums to zero. 




If you dig 
equations 
then you'll 



enjoy this article. If 
understanding how 
mathematical tools 
affect second-order 
systems is a source 
of much headscratch- 
ing for you, dig in and 
enjoy the insight 
Dave offers here. 
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a) 




b) 



1.5 2 2.5 
Time (s) 

Figure 1a— The aperiodic response of an overdamped, second- 
order system consists of the difference between two exponential 
decays, b— The periodic response of an underdamped, second- 
order system consists ofasinewave with an exponentially 
decaying amplitude. 
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solution can be written out 
like this: 



s 2 = -a- ya 2 -coo 2 
i(*)-AieV+ A#V 



A 1 and A 2 , 



The solution probably will be an expo- 
nential function. Let's try: 

i|t) = Ae s ' 

The first and second derivatives are: 



dt 



= Ase M 



dt 2 



= AsV 



Making the substitutions, the differen- 
tial equation becomes: 

^Ae" + RAse M +LAs 2 e sl =0 



Ae"(^ + Rs + Ls 2 )=0 



Solving this requires that the third 
factor in parentheses be set to zero. 
Because this is a quadratic equation, 
there are two solutions: 



s - — V 



-4L 



S = "2L + 



i 

LC 



By defining a couple of new constants, 
you can write a simpler form: 

°-2 R L and(0 °%LC 

a is the damping coefficient, and ca is 
the resonant frequency. In general, to 
is used to denote a radian frequency, 
which is related to a frequency / in 

www.circuitceitar.com 



hertz by a factor of 2k. I'm going to 
stick with radian frequency to keep 
the formulae simple; if you want to 
convert to hertz, just use the relation: 



2n 



You'll also find it useful to define the 
damping ratio: 

a. = r 

Look at the argument to the square 
root: If this quantity is negative, then 
the square root is an imaginary number 
and the solution s is going to be a 
conjugate pair of complex numbers. In 
real terms, this means that the decay is 
periodic. If it's exactly zero, then the 
circuit is said to be critically damped. 
And if the quantity is positive, there 
will be two real values of s in the 
solution and the decay will be aperi- 
odic. These three cases correspond to 
the quantity C, being less than, ex- 
actly, and greater than 1, respectively. 

It can be useful to 
solve this expression 
for how R relates to L 
and C (see initial qua- 
dratic solution): 



• if R 2 = 4L/C, the 
circuit is critically 
damped 

• if R 2 > 4L/C, the 
decay is aperiodic 

• if R 2 < 4L/C, the 
decay is periodic 

The full form of the 
CIRCUIT CELLAR* 



are picked to 
satisfy whatever initial 
conditions are set up for 
the circuit. If s 2 and s 2 axe 
two different real values, 
the overall function is the 
sum of two exponential 
equations with two differ- 
ent time constants. 
In general, you'll find 
that A 1 and A 2 have equal values but 
opposite signs, so you end up with a 
function that is zero at t = and at t = 
infinity (see Figure la). 

If s 1 and s 2 are a complex conjugate 
pair, you can define a different fre- 
quency, the natural frequency: 

o) N =v/oV-a 2 

And, this allows you to write the solu- 
tion as: 

i(t) = A l e(- a *»"N)' + A 2 e(-"-"»N)' 

Euler's Identity says: 

E jKl = cos Kt + j sin Kt 

Remembering that A, = -A 2 , you 
find that the cosine terms cancel and 
you're left with an expression that 
describes a decaying sinewave: 

i(t) = 2A,e^'j sin (coNt) 




-+- o 



Figure 2a— A sinewave can be characterized by its amplitude and phase angle, or 
equivalent^, by the magnitudes of its in-phase (sine) and quadrature (cosine) 
components. This vector can represent any sinewave quantity at a particular 
frequency. b-Current is in phase with the applied voltage in a resistor, but it lags 
by 90° in a coil and leads by 90° in a capacitor. 




Figure 3a— Here's the RLC circuit with a current source added in series. 
To find the total voltage, the component voltages must be added 
vectorially by placing the individual vectors end-to-end. b—At frequencies below 
resonance, the capacitive component is larger than the inductive component, 
resulting in an overall negative phase shift 0. c—At frequencies above reso- 
nance, the overall phase shift is positive. 

Figure lb shows an example of this 
type of response. 

The Mathcad sheet for Figure 1 lets 
you plug in values for R, L, and C to 
see the effect on a, co , co N , and £, 

DRIVEN RESPONSE 

This information provides insight 
into the basic mechanisms of second- 
order systems, but you're probably 
more interested in the response of the 
system to periodic driving forces. You 



did a simulation last 
month, now it's time 
to derive the math 
that describes this 
aspect of the system. 

This requires a 
different set of math 
tools from those used 
so far. Here, phasors 
and complex imped- 
ances rule the day, 
but first let me get a 
few basic concepts 
out of the way. 

Remember that 
the fast Fourier trans- 
form shows that any periodic forcing 
function can be expressed as a sum of 
sines and cosines, or equivalently, as 
a sum of sines whose relative phase 
can be varied. Because you're dealing 
with strictly linear circuits here, you 
can use the principle of superposition 
to analyze the circuit response one 
sinewave component at a time. And 
you can derive the overall response by 
summing the components at the end. 
Also, because of the linearity, the 



forced response must be a sinewave at 
the same frequency, but with a differ- 
ent amplitude and/or phase. 

Together these facts mean that if 
you understand the response of a cir- 
cuit or system to sinewaves of various 
frequencies, you always can get the 
response to arbitrary periodic forcing 
functions. Any sinewave is character- 
ized by just two parameters, its am- 
plitude and initial phase. Together, 
these parameters can be thought of as 
a vector in phase space, also known as 
a phasor (see Figure 2a). 

Amplitude can be in units of any 
physical quantity. This can let you 
show the relationship between volt- 
age and current for your components 
(see Figure 2b). In the resistor, the 
current always is in phase with the 
applied voltage; or, if you apply a 
current, the voltage is always in phase 
with that. The relationship between 
voltage and frequency is independent 
of the frequency. 

In the coil, the current always lags 
the voltage by 90°, or equivalently, 
the voltage leads the current by 90°. 
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Figure 4a— If the current source is replaced by a voltage source, the 
component relationships remain the same, but the diagrams must be 
rotated to put the overall voltage on the 0° axis, b— Below resonance, the current 
has a positive phase shift relative to the applied voltage. c—At a lower fre- 
quency, the current is smaller and the phase shift 9 is greater. 

However, their magnitudes are related 
by a function that depends on fre- 
quency. If the voltage is held con- 
stant, then the current decreases with 
increasing frequency, or if the current 
is held constant, the voltage increases 
with increasing frequency. Because 
impedance is defined to be the ratio of 
voltage to current, the last statement 
is the same as saying the impedance 
of the coil increases with frequency. 

Everything is reversed for the ca- 
pacitor. Thus, the voltage always lags 
the current by 90°, or equivalently, 
the current leads the voltage by 90°. 
Again though, their magnitudes are 
related by a function that depends on 
frequency. However, now, if the volt- 
age is held constant, the current esca- 
lates with increasing frequency, or if 
the current is held constant, the volt- 
age (and impedance) decreases with 
increasing frequency. 



Look at Figure 3a. 
Because the current 
is the same every- 
where in a series 
circuit, each of the 
components has a 
voltage across it 
determined by the 
effect of the current 
on that component. 
The voltage across 
the current source is 
the sum of the three 
component voltages. 
Figure 3b shows 
how they add up. Because the voltages 
are at different phase angles, use vec- 
tor addition to add them. The fre- 
quency of the current source is less 
than the resonant frequency of the 
circuit. The capacitor's voltage is 
greater than the coil's voltage, so the 
overall sum has a phase angle that's 
negative with respect to the current 
source (which is arbitrarily set at 0° 
for simplicity). 

As the frequency is decreased fur- 
ther, the C component continues to 
grow while the L component shrinks, 
causing the total vector to lengthen 
and the phase angle to increase. In the 
limit, the vector is infinitely long and 
parallel to the -90° axis. This corre- 
sponds to a DC current source charg- 
ing the capacitor indefinitely. 

Figure 3c shows the case where the 
source frequency is greater than the 
resonant frequency. Now, the coil's 
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Figure 5a — The sum of (he R and L voltage components pmduces a high-pass filter function, b— The output voltage 
and phase angle are given by the vector showing the partial sum of the R and L components. o—When the magni- 
tude and angle oftheR + L vector is plotted against frequency, you see the high-pass function. 
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Figure 6a— Swapping the R 
and L with the C produces a 
low-pass filter function. b— The 
magnitude of the C vector is 
the complement of the high- 
pass response, and the phase 
angles have the opposite 
values. 
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voltage component is greater than the 
capacitor's component (i.e., the coil's 
impedance is greater than the 
capacitor's, because the current is the 
same in both) and the phase angle of 
the sum is positive. If the frequency 
continues to increase, both the phase 
angle and amplitude increase, as well. 

As you might have guessed, the 
resonant case is the frequency at 
which the voltage (and impedance) of 
the coil and capacitor cancel each 
other, leaving only the resistive com- 
ponent in the sum. At this point, the 
total voltage is at a minimum and the 
phase angle is 0°, which means that the 
voltage is in phase with the current. 

What happens if a voltage source is 
substituted for the current source (see 
Figure 4a)? Well, you have the same 
sum of voltages as before, but now the 
total voltage is constrained to match 
that of the voltage source and the 
current phase angle and magnitude are 
free to vary. 

Figure 4b shows what happens 
when Figure 3b is rotated so that the 
sum voltage lies along the 0° axis, 
scaled so that the length of this vector 
corresponds to the applied voltage. 
The current is still the same every- 
where in the circuit, and because it is 
in phase with and proportional to the 
voltage across the resistor, just look at 
the resistive component of the voltage 
diagram to see what it is. 

The Mathcad sheet for Figures 3 
and 4 lets you to plug in various val- 
ues for R, L, C, and co to see the effect 
on the various component vectors. 



SECOND-ORDER FILTERS 

Second-order systems can be used as 
filters of various types. As you saw, the 
series RLC circuit has a notch in its 
voltage response when driven with a 
current source and a peak in its current 
response when driven with a voltage 
source. This gives basic band-stop and 
band-pass capabilities, but what about 
high-pass and low-pass functions? 
Start by redrawing the circuit of Fig- 
ure 4a (see Figure 5a), and take a look 
at the sum of the resistive and induc- 
tive components of the response. This 
is simply a matter of drawing an addi- 
tional vector in the phase sketch, as 
shown in Figure 5b, which shows the 
R + L component for several different 
frequencies. 

The curved, dotted line shows the 
locus of points reached by this vector 
for a wide range of frequencies. Note 
that w corresponds to point C on this 
curve, but this isn't where the V OUT 
vector is the longest, which actually 
occurs at point D. Point D corre- 
sponds to the co N calculated earlier. 

Figure 5c shows the length (magni- 
tude) and direction (phase angle) of 
the R + L vector plotted as a function 
of frequency. And sure enough, it 
looks like a high-pass filter. 

Go back to the physical mass- 
spring-damper system, where the 
spring corresponds to the capacitor, 
the mass to the coil, and the damper 
to the resistor. Recall that current in 
the electrical system corresponds to 
the motion of the mass (relative to 
the anchor point) in the physical sys- 



tern and that voltage corresponds to 
physical force. The AC voltage source 
in the electrical system is equivalent 
to a cyclical mechanical force applied 
to the anchor point of the mechanical 
system (which is free to move). 

If the cyclical frequency is high, 
the mass tends to sit still because of 
its inertia. The output voltage of the 
electrical system is equivalent to the 
total mechanical force across the 
damper and the net force on the mass. 
The damper resists high-speed mo- 
tion, so it couples most of the input 
force directly to the mass, and the 
anchor point moves very little. This 
means that there is insignificant force 
across the spring, which requires a 
strong physical motion of one end or 
the other to develop significant force. 
Any motion of the mass and anchor 
point will be in phase with each other 
and with the applied force. 

If the cyclical frequency is low, the 
mass moves in response to the forces 
on it, and because the damper provides 
little resistance to this low-velocity 
motion, most of this force is coupled 
by the spring. The anchor point moves 
as well, and the magnitude of the mo- 
tion of the mass relative to the anchor 
point (corresponding to current) is 
small. The damper plus mass force (the 
analog of the electrical output voltage) 
is negligible. 

The low-pass function is con- 
structed by realizing that whatever 
input voltage doesn't appear across 
the R-L combination must be across 
the C. So, if the R-L combination is 
swapped with the C (see Figure 6a), 
you get the inverse response shown in 
Figure 6b. You can do the same sort of 
vector analysis for the high-pass filter, 
but then the output vector is the same 
as the C component alone. 

UNTIL NEXT TIME 

Now you've seen some of the math 
that describes the transient behavior 
and steady state of second-order sys- 
tems. If you'd like to learn more, read 
the book listed in the Resources. It 
covers the topic with more detail and 
rigor, but in a very readable way. 

Next month, I'll talk about active 
second-order systems and then move 
onto the topic of servomechanisms. B 
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